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It  has  been  a  long-standing  goal  to  create  magnetism  in  a  non-magnetic  material  by 
manipulating  its  structure  at  the  nanoscale.  Many  structural  defects  have  unpaired  spins;  an 
ordered  arrangement  of  these  can  create  a  magnetically  ordered  state.  In  this  article  we  predict 
theoretically  that  stepped  silicon  surfaces  stabilized  by  adsorbed  gold  achieve  this  state  by  self- 
assembly,  creating  chains  of  polarized  electron  spins  with  atomically  precise  structural  order. 
The  spins  are  localized  at  silicon  step  edges  having  the  form  of  graphitic  ribbons.  The  predicted 
magnetic  state  is  supported  by  recent  experimental  observations,  such  as  the  coexistence  of 
double-  and  triple-period  distortions  and  the  absence  of  edge  states  in  photoemission.  Ordered 
arrays  of  surface  spins  can  be  accessed  by  probes  with  single-spin  sensitivity,  such  as  spin- 
polarized  scanning  tunnelling  microscopy.  The  integration  of  structural  and  magnetic  order  is 
crucial  for  technologies  involving  spin-based  computation  and  storage  at  the  atomic  level. 
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Individual  atoms  with  an  odd  number  of  electrons  have  a  mag¬ 
netic  moment  from  the  spin  of  the  unpaired  electron.  Materials 
consisting  of  elements  with  an  even  number  of  electrons — such 
as  carbon,  silicon  and  oxygen — can  also  reveal  unpaired  elec¬ 
trons  when  covalent  bonds  are  broken.  In  defective  or  disordered 
group  IV  solids,  the  single  electrons  that  occupy  defect  levels  can 
be  detected  in  magnetic  resonance  experiments1,  but  the  typically 
low  concentrations  of  these  defects  normally  preclude  a  magneti¬ 
cally  ordered  state.  Recent  intriguing  reports  of  ferromagnetism  in 
insulating  oxides  have  led  to  the  idea  of  ‘d°  magnetism’,  in  which 
magnetic  order  is  argued  to  arise  purely  from  spin-polarized  native 
defects  such  as  oxygen  or  cation  vacancies2.  These  phenomena 
remain  controversial  and  their  theoretical  understanding  is  incom¬ 
plete3'4.  Nevertheless,  the  pursuit  of  ferromagnetic  materials  without 
transition  metals  or  rare  earths  is  part  of  a  broad  effort  to  identify 
and  develop  new  systems  in  which  spin-polarized  electrons  can  be 
used  to  store  and  process  information. 

Bonds  are  also  broken  at  the  surfaces  and  interfaces  of  solids,  but 
covalent  materials  usually  reconstruct  to  eliminate  dangling  bonds  or 
doubly  occupy  them  with  electrons  of  opposite  spin.  Surfaces  none¬ 
theless  provide  natural  templates  for  creating  ordered  arrays  of  spins 
that  can  be  individually  addressed  by  local  probes.  For  example,  step 
edges  on  vicinal  Si(lll)  can  be  atomically  smooth  for  distances  as 
great  as  20,000  lattice  sites5.  These  step  edges  have  already  been  pro¬ 
posed  as  templates  for  creating  long  chains  of  nuclear  spins,  consist¬ 
ing  of  29Si  atoms,  as  the  basis  for  an  all-silicon  quantum  computer6. 

Here,  we  describe  a  very  different  alternative  to  fabricating  chains 
of  spins  on  templated  surfaces.  The  concept  rests  on  our  prediction 
that  a  family  of  well-known  vicinal  silicon  surfaces  are  themselves 
intrinsically  magnetic.  Thus,  these  surfaces  create — entirely  by  self- 
assembly — long  chains  comprising  thousands  of  electron  spins 
with  nearly  perfect  structural  order.  Such  ordered  arrays  of  spins 
in  silicon  systems  having  no  magnetic  elements  are  quite  unexpec¬ 
ted.  Indeed,  none  of  these  systems  has  yet  been  investigated  for 
magnetism  experimentally,  but  we  show  below  that  existing  data 
already  support  it  indirectly. 

Results 

Spin  polarization  at  graphitic  steps  in  silicon.  Our  investigation 
focuses  on  Si(553)-Au  and  Si(557)-Au  surfaces,  shown  in  Figure  1. 


These  are  stepped  surfaces  stabilized  by  a  fraction  of  a  monolayer 
of  Au,  and  are  members  of  a  family  of  Au-induced  vicinal  silicon 
surfaces  with  orientations  close  to  (111).  The  Si(553)-Au  surface 
is  the  most  perfect  of  the  Au-induced  chain  structures  on  vicinal 
Si(lll):  its  structural  perfection  is  limited  only  by  defects  created 
by  adsorption  of  residual  gas.  Even  these  will  be  frozen  out  at  the 
temperatures  normally  used  (below  5  K)  to  probe  individual  spins  at 
surfaces7.  Crystallographically,  Si(553)  and  (557)  can  be  considered 
miscut  from  the  flat  Si(  111)  surface  by  angles  of  12.3°  and  9.4°  away 
from  and  towards,  respectively,  the  (001)  orientation.  Several  other 
Au- stabilized  vicinal  (111)  surfaces  contain  structural  motifs  in 
common  with  Si(553)-Au  and  Si(557)-Au8.  These  commonalities 
suggest  that  magnetism  in  silicon  is  likely  to  be  found  throughout 
this  family. 

The  critical  structural  motif  for  magnetism  in  these  systems  is  a 
nanoscale  honeycomb  strip  of  graphitic  silicon,  just  a  single  hexa¬ 
gon  wide,  which  forms  the  edge  of  the  step9.  This  honeycomb  chain 
was  first  observed  in  reconstructions  of  flat  Si(  111)  induced  by  sub¬ 
monolayer  coverages  of  alkali,  alkaline  earth,  and  some  rare-earth 
metals  as  well  as  Ag  and  Au10.  In  these  flat  surfaces,  the  orbitals  of 
the  outer  honeycomb  atoms  form  filled  bands,  either  by  accept¬ 
ing  electrons  from  the  adsorbates  or  by  forming  partially  covalent 
bonds.  At  a  step  edge,  a  new  possibility  arises  in  which  these  orbitals 
are  singly  occupied,  fully  spin  polarized  and  magnetically  ordered. 

The  two  surfaces  studied  here  have  step  edges  with  opposite 
orientation,  as  shown  in  Figure  1.  We  used  density-functional  the¬ 
ory  (DFT)  to  determine  the  equilibrium  structure  of  these  surfaces. 
The  results  are  shown  in  Figure  la,b  for  Si(553)-Au  and  in  Figure  lc 
for  Si(557)- Au.  The  same  structural  model  for  Si(553)- Au  was  found 
independently  in  ref.  11  on  the  basis  of  non-spin-polarized  DFT 
calculations.  Similarly,  the  model  of  Si(557)-Au  is  well  established 
from  earlier  non-spin-polarized  DFT  calculations8,12  and  from 
X-ray  diffraction13. 

By  including  spin  polarization  in  our  calculations,  we  find  that 
some  of  the  dangling  bonds  at  the  step  edges  of  both  surfaces 
become  fully  polarized  with  a  spin  moment  of  one  Bohr  magneton. 
On  Si(553)-Au,  every  third  Si  atom  at  the  step  edge  is  completely 
spin  polarized,  whereas  all  other  atoms  are  negligibly  polarized.  On 
Si(557)-Au  (Fig.  lc),  every  second  Si  atom  at  the  step  is  fully  polar¬ 
ized,  as  are  all  the  terrace  restatoms.  In  both  systems,  the  interaction 


Figure  1  |  Ground  state  structure  and  lowest  energy  spin  configuration  of  two  magnetic  silicon  surfaces,  (a,  b)  Si(553)-Au  in  its  antiferromagnetic 
ground  state.  Yellow  atoms  are  Au,  all  others  are  Si.  Each  terrace  contains  a  Au  double  row  and  a  graphitic  Si  honeycomb  chain  (green)  at  the  step  edge. 
Every  third  Si  (red,  blue)  at  the  step  has  a  spin  magnetic  moment  of  one  Bohr  magneton  (S  =  1/2,  arrows)  from  the  complete  polarization  of  the  electron 
occupying  the  dangling-bond  orbital.  The  sign  of  the  polarization  alternates  along  the  step.  The  periodicity  along  the  step  is  tripled  by  a  small  downward 
displacement  of  the  magnetic  atoms.  The  periodicity  along  the  two  Au  rows  is  doubled,  with  alternating  short  and  long  bonds  a,.  The  full  1x6  unit  cell 
is  outlined,  (c)  Si(557)-Au  in  its  antiferromagnetic  ground  state.  The  step  direction  of  Si(557)  is  opposite  to  that  of  Si(553)  as  shown  by  the  schematic 
outlines.  Nevertheless,  this  surface  is  also  magnetic:  every  second  Si  (red,  blue)  at  the  step  has  a  magnetic  moment  of  one  Bohr  magneton,  as  do  all  the  Si 
restatoms  (blue,  red)  on  the  terrace. 
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between  neighbouring  spins  along  the  step  is  antiferromagnetic.  In 
the  spin-polarized  antiferromagnetic  ground  state,  the  total  energy 
of  the  system  is  lower  than  in  the  non-spin-polarized  state.  For 
Si(553)-Au,  this  energy  gain  is  20  meV  per  spin  (see  Supplementary 
Discussion  and  Supplementary  Fig.  SI),  suggesting  that  spin  polari¬ 
zation  will  be  important  only  well  below  room  temperature. 

For  Si(553)-Au,  the  spin-polarized  ground  state  solves  two  out¬ 
standing  experimental  puzzles:  (1)  scanning  tunnelling  microscopy 
(STM)  of  Si(553)- Au  reveals  the  coexistence  of  two  distortions  with 
different  periods  at  low  temperatures.  We  show  that  these  compet¬ 
ing  periodicities  are  not  only  correctly  predicted  but,  in  addition, 
are  inextricably  linked  with  the  spin-polarized  ground  state.  Their 
experimental  observation  constitutes  indirect  but  compelling  evi¬ 
dence  supporting  the  possibility  of  magnetism  in  Si(553)-Au.  (2) 
The  non-spin-polarized  band  structure,  although  largely  in  agree¬ 
ment  with  experiment  for  the  Au  chains,  predicts  sharp  step-edge 
states  not  seen  in  photoemission11.  These  ‘missing  states’  can  now 
be  understood  as  being  washed  out  by  fluctuations  in  the  exchange 
splitting  of  the  spin-polarized  electronic  structure. 

For  Si(557)-Au,  the  spin-polarized  calculation  resolves  a  differ¬ 
ent  experimental  puzzle:  the  observation  of  a  metal-insulator  tran¬ 
sition  at  low  temperature.  An  insulating  ground  state,  with  a  small 
gap,  is  indeed  found  in  the  spin-polarized  ground  state,  whereas  the 
non- spin-polarized  band  structure  is  metallic.  However,  the  most 
intense  energy  bands  in  photoemission  already  agree  well  with  the 
calculated  non-spin-polarized  bands  arising  from  the  Au  chains12. 
This  agreement  remains  equally  good  in  the  spin-polarized  descrip¬ 
tion.  In  both  descriptions,  a  small  band  splitting14-17  originates 
from  the  spin-orbit  interaction  of  the  Au  atoms12,18.  In  this  ‘Rashba 
effect’ — which  is  unrelated  to  the  magnetism  we  predict  here — 
the  two  split  bands  have  a  spin  polarization  of  100%  at  any  given 
momentum  in  reciprocal  space,  but  this  vanishes  when  integrated 
over  all  momenta  and  leaves  no  net  spin  polarization  in  real  space. 

Predictions  for  Si(553)-Au.  The  Si(553)-Au  surface  was  first 
observed  in  ref.  19  and  its  unusual  one-dimensional  electronic 
properties  have  since  been  extensively  studied.  Two  experimental 
findings  are  especially  striking:  (1)  STM  images  acquired  at  room 
temperature  show  alternating  bright  and  dim  rows,  each  with  the 
periodicity  of  the  Si  surface  lattice  constant  a020-23.  At  temperatures 
below  50  K  these  rows  separately  develop  higher-order  periodic¬ 
ity,  with  a  tripled  period  (3a0)  for  the  bright  rows  and  a  doubled 
period  (2 a0)  for  the  dim  rows.  These  effects  have  been  attributed  to 
periodic  lattice  distortions21  and  to  charge-density  waves22,  but  their 
true  microscopic  origin  has  remained  elusive.  (2)  Angle-resolved 
photoemission  spectra  (ARPES)  obtained  at  160-220  K  reveal  three 
metallic  bands,  all  with  their  minima  at  the  zone  boundary  (in  the 
direction  parallel  to  the  chains)  ZB,,  of  the  lxl  Brillouin  zone19. 
The  two  lower  bands  are  slightly  split  by  a  small  momentum  split¬ 
ting  5kF  =  0.05  A-1  (ref.  18),  and  cross  the  Fermi  level  about  halfway 
to  the  zone  centre.  A  detailed  analysis  of  this  crossing  shows  that 
the  splitting  is  due  to  the  spin-orbit  interaction18.  The  upper  band 
crosses  the  Fermi  level  about  one-fourth  of  the  way  to  the  zone  cen¬ 
tre  without  a  detectable  splitting. 

All  these  observations  are  correctly  predicted  by  the  proposed 
model  of  Si(553)-Au.  Figure  2  shows  the  theoretically  simulated 
STM  image  for  the  model  in  Figure  la  in  its  antiferromagnetic 
ground  state.  The  bright  row  arises  from  Si  atoms  at  the  step  edge 
and  indeed  has  tripled  periodicity  3 a0.  The  highest  peaks  within 
this  row  are  from  the  spin -polarized  atoms.  This  appears  counter¬ 
intuitive,  as  these  atoms  are  actually  0.3  A  below  the  height  of  their 
non-polarized  neighbours.  But  each  spin-polarized  atom  has  one 
unoccupied  orbital  just  above  the  Fermi  level,  whereas  the  non¬ 
polarized  atoms  have  only  occupied  orbitals.  For  tunnelling  into 
empty  states,  the  topographic  peaks  are  therefore  indeed  higher  for 
spin-polarized  atoms. 
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Figure  2  |  Scanning  tunnelling  microscopy  of  Si(553)-Au.  (a)  Theoretical 
STM  topography  showing  the  tip  height  at  constant  current  for  tunnelling 
into  empty  surface  states  at  bias  voltage  +0.5V.  Inset:  electron  spin 
density  1 A  above  the  step  edge.  Red  and  blue  contours  represent  positive 
and  negative  polarization,  indicating  antiferromagnetic  alignment  of  the 
red  and  blue  Si  atoms  the  step  edge.  Contours  are  separated  by  a  factor 
of  two.  (b)  Line  scans  at  the  locations  marked  in  a,  showing  coexisting 
3a0  and  2o0  periodicities  along  the  Si  step  edge  and  Au  chain,  respectively. 


The  dim  second  row  visible  in  Figure  2a  arises  from  the  right- 
hand  leg  of  the  ladder  formed  by  two  rows  of  Au  atoms.  Its  doubled 
periodicity  2a0  is  due  to  the  dimerization  created  when  the  rungs 
of  this  ladder  rotate  with  alternating  signs  along  a  row.  A  similar 
dimerization  occurs  in  a  closely  related  system,  Si(lll)-(5x2)-Au, 
where  it  is  driven  by  electron  doping  from  Si  adatoms  and  the 
consequent  formation  of  a  surface  band  gap24.  For  Si(lll)-(5x2)-Au, 
the  dimerization  parameter  d=(a1-fl0)/a0,  where  a0  is  the  surface 
lattice  constant  and  a1  is  defined  in  Figure  lb,  has  the  value  d  =  0.14, 
much  larger  than  the  value  d  =  0.04  for  Si(553)-Au.  But  even  this 
small  dimerization  is  sufficient  to  produce  a  doubled  periodicity 
that  is  well  resolved  in  the  theoretical  STM  image. 

The  spin  polarization  of  the  electron  density,  shown  in  the  inset  of 
Figure  2a,  is  strongly  localized  at  every  third  Si  step-edge  atom.  Inte¬ 
grating  this  spin  density  shows  that  the  red  and  blue  step-edge  atoms 
each  have  a  spin  moment  of  one  Bohr  magneton,  equivalent  to  a  fully 
polarized  electron  with  spin  S  =  1/2.  The  direction  of  the  spin  alter- 
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nates  along  the  step  in  this  antiferromagnetic  ground  state.  Moments 
on  other  atoms  along  the  step  are  smaller  by  an  order  of  magnitude. 

The  connection  between  magnetism  and  the  periodicity  of  the 
rows  is  simple  and  direct.  If  the  spin  polarization  is  constrained  to 
be  zero,  the  DFT  equilibrium  geometry  of  the  model  changes:  both 
the  tripled  periodicity  along  the  Si  rows  and  the  doubled  periodic¬ 
ity  along  the  Au  rows  are  completely  eliminated.  In  this  sense,  the 
experimental  observation  of  triple  and  double  periodicities  consti¬ 
tutes  strong  evidence  for  a  spin-polarized  ground  state. 

Magnetism  also  provides  the  key  to  understanding  a  long-stand¬ 
ing  puzzle  presented  by  ARPES  data  for  Si(553)-Au.  The  orbitals  at 
the  step  edges  of  a  vicinal  surface  would  normally  give  rise  to  sharp 
electronic  states  near  the  Fermi  level,  as  is  indeed  the  case  for  previ¬ 
ously  proposed  models  of  Si(553)-Au25,26.  But  these  are  not  observed 
in  the  ARPES  data.  We  demonstrate  below  that  spin  polarization  at 
the  step  edge  provides  the  resolution  to  this  puzzle. 

Figure  3a  shows  the  theoretical  bands  for  Si(553)-Au  unfolded 
into  the  extended  Brillouin  zone,  where  they  can  be  easily  compared 
with  those  obtained  from  ARPES  spectra.  As  the  theoretical  ground 
state  is  antiferromagnetic,  one  need  only  analyse  the  bands  from 
a  single  spin  channel.  The  symbols  and  curves  identify,  using  the 
colours  from  Figure  1,  the  states  having  the  largest  projection  onto 
particular  atoms  at  the  surface.  States  arising  from  the  magnetic 
atoms  experience  an  exchange  splitting  of  about  0.5  eV.  As  a  result, 
for  spin-up  states,  the  red  step-edge  band  is  singly  occupied  and  the 
blue  step-edge  band  is  empty.  For  spin-down  states,  this  situation  is 
reversed,  as  seen  in  the  density  of  states  (DOS;  Figure  3b).  A  more 
accurate  treatment  of  electron  correlation  (beyond  DFT)  could  shift 
these  narrow  bands  somewhat,  but  the  large  exchange  splitting  sug¬ 
gests  that  qualitative  changes  are  unlikely. 

Although  these  magnetic  step-edge  states  are  responsible  for  the 
1x6  structural  and  magnetic  periodicity  of  Si(553)-Au,  the  remain¬ 
ing  features  of  the  band  structure  have  a  simple  lxl  non-magnetic 
character.  This  is  evident  from  comparing  the  1x6  bands  with  those 
in  Figure  3d  for  a  simpler  lxl  structure  identical  to  our  model  but 
constrained  to  unit  periodicity  a0  and  thus  having  no  spin  polariza¬ 
tion.  In  this  lxl  structure,  there  is  just  a  single  type  of  step  atom  and 
therefore  only  one  corresponding  band  (green).  The  remnant  of  this 
non-magnetic  band  is  also  visible  in  the  full  1x6  bands,  although 
there  its  spectral  weight  at  a  given  momentum  is  spread  out  by 
partial  hybridization  with  its  magnetic  counterparts. 


Both  the  1x6  and  lxl  structures  give  rise  to  several  metallic 
bands  (yellow)  with  their  minima  at  the  zone  boundary  ZBy.  These 
are  the  bands  detected  in  ARPES  (Fig.  3d,  inset).  In  the  1x6  case, 
there  are  four  bands  in  the  two  spin  channels,  but  the  spin-up  and 
spin-down  bands  are  nearly  degenerate.  In  the  non-magnetic  lxl 
case,  there  are  also  four  bands,  but  here  the  degeneracy  is  bro¬ 
ken  by  the  spin-orbit  interaction27,  which,  for  simplicity,  was  not 
included  in  the  1x6  calculation.  The  resulting  momentum  splitting 
5kF  =  0.04  A-1  of  the  lower  band  is  in  excellent  agreement  with  the 
experimental  splitting.  The  splitting  of  the  upper  band  is  less  than 
half  this  value  and  consequently  is  not  resolved  in  the  ARPES  data. 
The  predicted  Fermi  level  crossings  of  these  bands  are  in  qualitative 
agreement  with  the  measured  values:  the  lower  doublet  crosses  the 
Fermi  level  halfway  to  the  zone  centre,  and  the  upper  doublet  about 
one-fifth  of  the  way. 

One  feature  of  the  non-magnetic  lxl  band  structure  is  not 
observed  in  the  ARPES  data:  the  green  step-edge  band  that  just 
crosses  the  Fermi  level  in  Figure  3d.  Magnetism  provides  a  simple 
explanation  for  the  apparent  absence  of  this  band.  Above  a  tempera¬ 
ture  scale  set  by  the  energy  gained  from  spin  polarization,  thermal 
fluctuations  between  polarized  and  unpolarized  electron  states  will 
broaden  the  step-edge  bands.  The  magnitude  of  the  broadening  is 
roughly  the  exchange  splitting,  0.5  eV.  The  resulting  reduction  in 
intensity  makes  the  step-edge  bands  essentially  invisible  compared 
with  the  sharp  Au  bands,  which  are  not  affected  by  the  spin  polari¬ 
zation  at  the  step  edge.  In  this  sense  the  experimental  absence  of 
a  sharp  step-edge  band  constitutes  additional  evidence  for  a  spin- 
polarized  ground  state. 

Predictions  for  Si(557)-Au.  We  now  briefly  summarize  our  predic¬ 
tions  for  the  magnetic  properties  of  Si(557)-Au.  The  main  features  are 
similar  to  those  of  Si(553)-Au,  with  several  interesting  differences. 

(1)  Crystallographically,  the  steps  on  Si(557)  are  oriented  in  an 
opposite  manner  to  those  on  Si(553).  On  clean,  unreconstructed  sil¬ 
icon,  this  difference  would  lead  to  two  dangling  bonds  at  each  step- 
edge  atom,  in  contrast  to  the  single  dangling  bond  at  each  step-edge 
atom  of  Si(553).  But  the  graphitic  reconstruction  of  the  steps  on 
Si(557)-Au  leads,  just  as  for  Si(553)-Au,  to  a  single  dangling  bond 
per  step-edge  atom.  We  find  that  every  second  of  these  dangling 
bonds  is  fully  spin  polarized.  This  2 a0  periodicity  in  spin  polariza¬ 
tion  tracks  the  up-down  buckling  of  the  step-edge  atoms  previously 
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Figure  3  |  Electronic  band  structure  and  density  of  states  (DOS)  for  Si(553)-Au.  (a,  b)  Bands  and  DOS  for  the  ground-state  antiferromagnetic  structure 
of  Figure  la.  The  full  1x6  bands  are  shown  unfolded  and  parallel  to  the  chain  direction  (ZB||  is  the  boundary  of  the  1x1  Brillouin  zone)  along  the  line  halfway 
to  the  orthogonal  zone  boundary  ZB±.  Red  and  blue  diamonds  are  spin-polarized  states  from  red  and  blue  Si  step-edge  atoms  in  Figure  la.  Only  spin-up 
states  are  shown  here— antiferromagnetic  order  implies  the  spin-down  states  are  equivalent  but  with  red  and  blue  reversed,  as  seen  in  the  DOS.  Green 
diamonds  are  non-polarized  states  from  green  Si  step-edge  atoms.  Yellow  circles  are  bonding  (lower)  and  anti-bonding  (upper)  states  from  the  Au  chains. 
Curves  show  schematic  energy  bands  without  hybridization.  Projected  bulk  states  are  not  shown,  for  clarity,  (c,  d)  Bands  and  DOS  for  the  same  model 
when  it  is  constrained  to  1x1  periodicity,  for  which  the  ground  state  is  not  spin-polarized.  Spin-orbit  coupling,  included  in  the  bands  but  not  in  the  DOS, 
splits  each  of  the  two  Au  chain  states  but  not  the  Si  step-edge  state.  Inset:  experimental  angle-resolved  photoemission  spectra  from  Crain  et  al.19 
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found  in  non-spin-polarized  calculations12,28.  Only  the  down  atoms 
are  polarized. 

(2)  The  terraces  on  Si(557)-Au  are  wider  than  those  on  Si(553)- 
Au  and  contain  a  row  of  ‘restatoms’  not  found  on  Si(553)-Au.  We 
find  that  the  dangling  bond  of  each  restatom  is  singly  occupied  and 
fully  spin  polarized.  Because  the  restatoms  are  separated  by  2 a0 
along  the  step  direction,  the  number  of  spin-polarized  step-edge 
atoms  and  restatoms  is  the  same. 

(3)  The  dominant  interaction  between  neighbouring  spins  along 
the  step  edge  is  antiferromagnetic;  the  same  is  true  for  the  spins 
along  the  row  of  restatoms.  Adjacent  rows  of  spins  are  also  prefer¬ 
entially  aligned  antiparallel.  In  this  antiferromagnetic  ground  state, 
the  energy  of  the  system  is  reduced  by  45  meV  per  spin  relative  to 
the  non-spin-polarized  state. 

(4)  In  non-spin-polarized  calculations,  the  restatom  orbital 
is  doubly  occupied  and  the  orbital  of  the  down  step-edge  atom  is 
nearly  empty12,28.  In  the  spin-polarized  ground  state,  this  changes: 
both  the  restatom  and  down  step-edge  atoms  are  singly  occupied 
(and  polarized).  This  charge  transfer  reduces  the  magnitude  of  the 
step-edge  buckling  and  also  reduces  the  height  of  the  restatoms. 

(5)  There  is  considerable  experimental  controversy  regarding  the 
possibility  that  Si(557)-Au  undergoes  a  metal-insulator  transition 
at  low  temperature17,29-31.  The  origin  of  an  insulating  state  at  low 
temperature  is  not  clear.  The  theoretical  non-spin-polarized  band 
structure  is  metallic.  It  has  been  proposed  that  thermal  fluctuations 
at  the  step  edge  could  contribute  to  a  metal-insulator  transition12,28. 
An  alternative  explanation  is  found  in  the  spin-polarized  band 
structure,  which  is  insulating  with  a  small  gap  of  75meV  (which 
persists  when  spin-orbit  coupling  is  included). 

(6)  The  first  theoretical  studies  of  Si(557)- Au  included  spin-orbit 
coupling  but  not  spin  polarization12.  The  resulting  prediction  of 
momentum- split  Au  bands  was  in  good  agreement  with  experi¬ 
ments.  These  studies  also  predicted  Si  step-edge  and  restatom  bands 
not  seen  in  ARPES8,15,17.  In  Riikonen  et  al.2S  it  was  demonstrated 
theoretically  that  thermal  fluctuations  cause  these  bands  to  oscil¬ 
late  with  an  amplitude  of  approximately  0.2  eV  and  can  thus  blur 
these  nominally  sharp  states.  By  including  spin  polarization  in  the 
calculation,  another  scenario  emerges.  The  exchange  splitting  in 
Si(557)-Au  pushes  the  two  occupied  spin  bands  down  to  0.4-0. 5 eV 
below  the  Fermi  level,  and  the  two  unoccupied  spin  bands  up  to 
0.1-0.2eV  above  the  Fermi  level.  This  exchange  splitting,  approxi¬ 
mately  0.5 eV,  is  similar  to  that  we  found  for  Si(553)-Au.  Thermal 
fluctuations  between  polarized  and  unpolarized  states  will  similarly 
broaden  these  step-edge  and  restatom  bands,  rendering  them  invis¬ 
ible  to  ARPES  while  leaving  the  momentum- split  Au  bands  unaf¬ 
fected.  The  scale  of  the  fluctuations  in  both  scenarios  is  comparable, 
and  thus  it  is  likely  that  both  mechanisms  contribute  to  blurring  the 
energy  bands. 

(7)  In  contrast  to  our  findings  for  Si(553)-Au,  we  do  not  predict 
any  higher-order  periodicity  to  mark  the  onset  of  magnetism  at  low 
temperature:  the  periodicity  of  Si(557)-Au  remains  1x2  in  both  the 
non- spin-polarized  and  spin-polarized  states. 

Discussion 

We  anticipate  that  spin-polarized  states  in  silicon  are  not  limited  to 
Si(553)-Au  and  Si(557)-Au,  because  the  key  structural  element — a 
graphitic  step  edge — is  known  to  exist  on  other  stepped  silicon  sur¬ 
faces  as  well8.  Direct  experimental  tests  of  our  predictions  using,  for 
example,  spin-polarized  STM32,33  or  spin-polarized  photoemission34 
will  be  an  important  next  step,  and  will  provide  important  informa¬ 
tion  about  the  nature  of  the  magnetic  ordering. 

To  set  the  stage  for  such  tests,  it  is  helpful  to  estimate  the  mag¬ 
nitude  of  the  spin  interactions  within  a  simple  nearest-neighbour 
Heisenberg  Hamiltonian.  From  DFT  total-energy  calculations  of  dif¬ 
ferent  spin  configurations  on  Si(553)-Au,  we  find  antiferromagnetic 
coupling  (/||=15meV)  along  the  steps  and  weaker,  ferromagnetic 


coupling  (/±  =  -  0.3  meV)  across  the  steps  (these  values  may  be  dif¬ 
ferent  in  methods  that  treat  electron  correlation  explicitly;  see  the 
discussion  in  Volnianska  et  al.3S).  Intriguingly,  we  find  that  the  mag¬ 
nitude,  and  even  the  sign,  of  these  couplings  can  be  changed  by  dop¬ 
ing  electrons  or  holes  into  the  Si(553)-Au  surface  states  (see  Supple¬ 
mentary  Methods  and  Supplementary  Fig.  S2).  It  is  well  established 
that  surface  states  on  the  closely  related  Si(lll)-Au  surface  can 
be  electron-doped  by  adsorbates  (for  example,  silicon  adatoms) 
on  the  surface,  and  that  by  varying  this  adsorbate  population  one 
can  perform  band- structure  engineering  with  extraordinary  preci¬ 
sion36.  The  possibility  of  tuning  surface  magnetism  using  adsorbate 
dopants,  or  surface  chemistry  more  generally,  suggests  a  fascinating 
new  research  direction. 

Linear  chains  of  spin-polarized  atoms  provide  atomically  perfect 
templates  for  the  ultimate  memory  and  logic,  in  which  a  single  spin 
represents  a  bit.  Here,  the  spins  are  locked  into  a  self-assembled 
rigid  lattice  on  a  substrate  that  is  compatible  with  silicon  technol¬ 
ogy.  One  potential  application  is  the  spin  shift  register  recently  pro¬ 
posed  theoretically  by  Mahan37.  This  device  requires  a  one-dimen¬ 
sional  chain  of  identical  atoms,  each  with  spin  S  =  1/2.  When  an 
additional  electron  is  conducted  along  the  chain,  each  spin  state  is 
shifted  by  one  atom — the  spin  analogue  of  a  standard  shift  register 
memory  device.  The  arrangement  of  spins  is  arbitrary  and  there  is 
no  requirement  of  long-range  order.  However,  two  other  important 
criteria  must  be  satisfied:  in  the  ground  state  each  atom  must  have 
one  electron,  and  the  atoms  must  have  correlated  electronic  states. 
The  first  criterion  is  indeed  satisfied  by  the  spin-polarized  atoms  on 
Si(553)-Au  and  Si(557)-Au.  In  previous  research  we  showed  that 
the  electronic  states  of  singly  occupied  silicon  dangling  bonds  are 
indeed  strongly  correlated38. 

Another  potential  application  is  the  storage  of  information  in 
single  magnetic  atoms39.  Towards  this  end,  Hirjibehedin  et  al.40 
recently  used  spin  excitation  spectroscopy  to  demonstrate  that  the 
orientation  of  an  individual  spin  can  be  measured  with  an  STM. 
For  information  storage,  the  stability  of  this  orientation,  for  example 
from  coupling  of  the  spin  to  the  substrate,  is  critical.  This  coupling  is 
most  likely  to  be  substantial  when  the  substrate  is  anisotropic  at  the 
atomic  scale.  In  the  experiments  by  Hirjibehedin,  this  anisotropy 
was  engineered  by  creating  Cu2N  islands  that  break  the  four-fold 
symmetry  of  the  Cu(001)  substrate.  In  Si(553)-Au  and  related  vici¬ 
nal  silicon  surfaces,  this  anisotropy  arises  naturally  at  the  step  edge. 
The  investigation  of  spin-lattice  coupling  in  these  magnetic  sili¬ 
con  surfaces  is  an  important  next  step  towards  the  ultimate  goal  of 
spintronics  with  single  spins. 

Methods 

Density-functional  theory  calculations.  First-principles  total-energy  calculations 
were  used  to  determine  equilibrium  geometries  and  relative  energies  of  the  models 
in  Figure  1  and  their  variants.  The  calculations  were  performed  in  a  slab  geometry 
with  six  or  more  layers  of  Si  along  with  the  reconstructed  top  surface  layer  and  a 
vacuum  region  of  at  least  10  A.  All  atomic  positions  were  relaxed,  except  the  bot¬ 
tom  Si  layer  and  its  passivating  hydrogen  layer,  until  the  largest  force  component 
on  every  atom  was  below  0.02  eV/A.  Total  energies  and  forces  were  calculated 
within  the  generalized-gradient  approximation  of  Perdew,  Burke  and  Ernzerhof41 
to  DFT  using  projector- augmented  wave  potentials,  as  implemented  in  VASP42,43. 
The  plane-wave  cutoff  for  all  calculations  was  250  eV. 

Brillouin-zone  sampling.  The  sampling  of  the  surface  Brillouin  zone  was  chosen 
according  to  the  size  of  the  surface  unit  cell  and  the  relevant  precision  require¬ 
ments.  For  the  1x6  reconstruction  of  Si(553)-Au  shown  in  Figure  la,b,  we  calcu¬ 
lated  the  equilibrium  geometries  and  total  energies  of  non-magnetic,  ferromagnetic 
and  antiferromagnetic  configurations  using  6x6  zone  sampling,  with  convergence 
checks  using  10x10  sampling.  For  the  1x2  reconstruction  of  Si(557)-Au  shown 
in  Figure  lc,  we  calculated  the  equilibrium  geometries  and  total  energies  of  non¬ 
magnetic,  ferromagnetic  and  antiferromagnetic  configurations  using  2x8  zone 
sampling,  with  convergence  checks  using  2x16  sampling. 

Exchange  coupling.  Exchange  coupling  constants  /y  and  J±  for  Si(553)-Au  were 
extracted  from  DFT  total  energies.  Supercells  containing  two  spins  on  each  of  two 
adjacent  steps  were  constructed  with  the  four  possible  unique  spin  configurations. 
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The  DFT  total  energies  were  fit  to  a  nearest-neighbour  classical  Heisenberg 
Hamiltonian.  The  numerical  reliability  of  the  coupling  constants  was  analysed  with 
respect  to  Brillouin-zone  sampling,  plane-wave  cutoff  and  relativistic  treatment, 
and  the  dependence  of  the  coupling  constants  on  charge  doping  of  the  surface 
states  was  investigated  (see  Supplementary  Methods  and  Supplementary  Tables). 

Scanning-tunnelling  microscopy  imagery.  The  simulated  STM  image  in  Figure 
2  was  calculated  using  the  method  of  Tersoff  and  Hamann44.  For  this  empty-state 
image,  we  integrated  the  local  density  of  states  (LDOS)  over  an  energy  window 
of  unoccupied  states  up  to  0.5  eV  above  the  Fermi  level.  The  simulated  STM  topo¬ 
graphy  under  constant- current  conditions  was  obtained  by  plotting  the  height  at 
which  the  integrated  LDOS  is  constant. 
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